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ABSTRACT 

We incorporate a contribution to reionization from X-rays within analytic and 
semi-numerical simulations of the 21-cm signal arising from neutral hydrogen during 
the epoch of reionization. The relatively long X-ray mean free path (MFP) means 
that ionizations due to X-rays are not subject to the same density bias as UV ion- 
izations, resulting in a substantive modification to the statistics of the 21-cm signal. 
We explore the impact that X-ray ionizations have on the power spectrum (PS) of 
21-cm fluctuations by varying both the average X-ray MFP and the fractional contri- 
bution of X-rays to reionization. In general, prior to the epoch when the intergalactic 
medium is dominated by ionized regions (Hii regions). X-ray- induced ionization en- 
hances fluctuations on spatial scales smaller than the X-ray MFP, provided that X-ray 
heating does not strongly supress galaxy formation. Conversely, at later times when 
Hii regions dominate, small-scale fluctuations in the 21-cm signal are suppressed by 
X-ray ionization. Our modelling also shows that the modification of the 21-cm signal 
due to the presence of X-rays is sensitive to the relative scales of the X-ray MFP, 
and the characteristic size of Hii regions. We therefore find that X-rays imprint an 
epoch and scale-dependent signature on the 21-cm PS, whose prominence depends 
on fractional X-ray contribution. The degree of X-ray heating of the IGM also de- 
termines the extent to which these features can be discerned. We further show that 
the presence of X-rays smoothes out the shoulder-like signature of H ii regions in the 
21-cni PS. For example, a 10% contribution to reionization from X-rays translates to 
a 20-30% modulation in the 21-cm PS across the scale of Hii regions. We show that 
the MWA will have sufficient sensitivity to detect this modification of the PS, so long 
as the X-ray photon MFP falls within the range of scales over which the array is most 
sensitive (~ O.lMpc"^). In cases in which this MFP takes a much smaller value, an 
array with larger collecting area would be required. As a result, an X-ray contribution 
to reionization has the potential to substantially complicate analysis of the 21-cm PS. 
On the other hand, a combination of precision measurements and modelling of the 
21-cm PS promises to provide an avenue for investigating the role and contribution of 
X-rays during reionization. 

Key words: cosmology: diffuse radiation, theory - galaxies: high redshift, intergalac- 
tic medium 



1 INTRODUCTION 

The process of hydrogen reionization started with ionized 
(Hii) regions around the first galaxies, which later grew to 
surround groups of galaxies. Reionization completed once 
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these H ii regions overlapped and occupied most of the vol- 
ume between galaxies. The inhomogeneous nature of the 
reionization process suggests that the spatial dependence 
of the statistics of redshifted 21-cm fluctuations will pro- 
vide the most powerfu l probe of the reionization epoch 
jPurlanetto et al.ll2006l ). With this as motivation, several 
experiments are currently under development that aim to 
detect the 21-cm signal during reionization, including the 
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Low Frequency ArrajQ (LOFAR) , the Murchison Widefield 
Arraj0 (MWA) and the Precision Array to Probe Epoch 
of ReionizatioiJl (PAPER), and more ambitious designs 
are being planned such as the Square Kilometre ArrajQ 
(SKA). Detection of the redshifted 21-cm signal will not only 
probe the astrophysics of reionization, but also the mat- 
ter power spectrum (PS) during the epoch of reionization 
l|McQuinn et allbood : [Bowman et al.ll2006al ). 

In anticipation of forthcoming 21-cm observations, 
much recent theoretical attention has focused on the 
prospects of measuring the PS of 21-cm emis si on (e.g. 



Zaldarriaga et al. 2004 



. .Furlanetto et al.1 12004 iMoralesI 
2005; Bowm an et al.1 l2006bt ). In particular, the ionization 
structure of the intergalactic medium (IGM) owing to UV 
emission associated with st ar formation has been s tudied in 
detai l using analytic (e.g. iFurlanetto et al.l |2004|: iBarkanal 
20071) and nume rical simulations (e.g. McQuinn et al. 20061 : 



Iliev et al]|2008l ) . These studies describe a scenario in which 



very large H II regions form around clustered sources within 
overdense regions of the IGM. The formation of these Hii 
regions has a significant effect on the shape of the 21-cm 
PS, moving power from small to large scales and thus leav- 
ing a shoulder-shaped feature on the PS at the char acter- 
istic scale of the Hll regions (|Furlanetto et al.l 120041 '). The 
detailed morphology of the ionization structure, which is 
encoded in the PS, will yield information about the ion- 
izing sources, as well as the structure of absorbers in the 
IGM on small scales (|McQuinn et al.ll200d ). More recently, 
semi-numerical models have also been developed to study 
the formation of HlI regions during the epoch of reioniza- 
tion (EoR). These have again focu sed on the role of UV 
emission from star- forming regio ns (Mesi nger fc Furlanetto! 
l2007l ) and quasars (|Geil fc Wvitha,2008i ). 

In addition to UV photons from the first galax- 
ies and quasars, it has been suggested that a back- 
ground of X-ray photons at very high redshift may 
be imp ortant in explaining the large Thomson optical 
depth (iDunklev et al.1 120081) of t he IGM measured by 
iyM_ AP (|Ricotti fc Ostrikerl |2004 IPritchard fc Furlanettol 
l2007l ). Indeed, several authors have interpreted this ob- 
servation by proposing an initial phase of preheating and 
partial ionization of the IGM by X-rays, resulting pri- 
maril y from black hole acc r etion (jShull fc van Steenberd 



19851: iMadau fc Efstathioul 1 19991: iHaiman et al.1 l2000l : 



Venkatesan et all 1200 ll : iRicotti fc Ostrikeill2004h . An X-ray 



background could also have originated from X-ray binaries 
and supernova remnants, through accretion in the former 
case and inverse Comp ton scattering in the latter (see, e.g. 
IFurlanetto et al]|2006l . for a review of these processes and 
the global evolution of the IGM). 

The process by which this early X-ray background ion- 
izes and heats the IGM begins with the photoionization of 
hydrogen and helium. The resulting high energy free elec- 
trons then deposit their energy through coUisional ioniza- 
tion and excitation of hydrogen and helium, and electron- 
electron scattering with cooler free electrons. Thus X-rays 
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ionize hydrogen both directly and through secondary ioniza- 
tions by photoelectrons from ionized helium, with the latter 
dominating so that many hydrogen atoms can be ionized by 
a single X-ray photon. As reionization proceeds, however, 
the proportion of each X-ray's energy that is deposited into 
the IGM as heat, increases. In particular, for ionization frac- 
tions > 1 0%, ionization by second ary electrons becomes in- 
efficient l|Ricotti fc Ostrikerl |2004i ). as the electron-electron 
scattering cross-section becomes large. As a result. X-ray 
ionization is self regulating at the level of ~ 10 — 20%. Fur- 
th er d e tails o f the X-ray emission mechanisms are discussed 
in lOhl (|200ll ). In addition, observational limits on the pos- 
sible contribution of X-rays to reioniz ation are provided by 
the unresolved soft X-ray background. iDiikstra et al.l (|2004l ) 
show that the hard X-ray background produced by a popu- 
lation of black holes at high redshift would be observed as 
a present-day soft X-ray background. They find that accret- 
ing black holes with a hard spectrum could not reionize the 
Universe to greater than the 50% level without saturating 
the unresolved component of the soft X-ray background. 

IPritchard fc Furlanettol l|2007f) study the eff ect of 

the inhomogeneous preheating of the IGM (|Ohl I2OOII : 
I Venkatesan et al.ll200ll : lMadau et al.ll2004l ) that resuhs from 
fluctuations in the X-ray background, and calculate its im- 
pact on 21-cm fluctuations at 2 ^ 10. The effect of this X-ray 
heatin g has also been ca lculated within numerical simula- 
tions ( Santos et al.|[2007l ). Our work is distinct from these 
studies, which do not include an X-ray component of reion- 
ization. In our modelling we neglect temperature fluctua- 
tions, and instead concentrate on the possible effect of X- 
ray ionization late in the reionization era when H II regions 
dominate the ionization structure, and the IGM is already 
heated above the cosmic microwave background (CMB) tem- 
perature. 

The Hll region dominated ionization structure of the 
IGM seen in numerical simulations forms as a result of the 
biased clustering of sources, combined with the short mean 
free path (MFP) of UV photons. In addition to the rela- 
tively large amount of energy deposited into the IGM, X-rays 
differ from UV photons in their effect on the IGM due to 
their comparatively long MFP. This long MFP means that 
X-rays provide a weakly fluctuating radiation background, 
rather than the highly biased and strongly nonlinear ionizing 
structure thought to arise from UV photons. The ionization 
structure of the IGM would therefore be modified if X-rays 
contributed significantly to reionization, because the long 
X-ray MFPs imply that their ionizing effect is not confined 
to th e vic inity of large scale overdensities in the matter field 
(Oh 20Q]I)- 

For this reason, the inclusion of an X-ray compo- 
nent of the emission responsible for the reionization of the 
IGM could modify fiuctuations in the 21-cm PS by reducing 
the influence of strongly biased ionized regions. 

In this paper we extend existing analytic and semi- 
numerical simulations of reionization to include an X-ray 
contribution to the ionizing flux. The structure of our pa- 
per is as follows: we begin by summarising estimates of the 
relative X-ray background (§[2| arising from stellar and ac- 
creting stellar mass black hole sources. §[3] describes our an- 
alytic model of the 21-cm signal arising from the EoR. In 
§|4]we describe how this model is extended to include a flux 
contribution from X-rays, and present results for the effect 
of X-ray ionization on the fluctuation statistics. §[5]summa- 
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rizes the semi-numerical scheme used to generate ionization 
maps of the IGM arising from both UV and X-ray emission. 
Our semi-numerical results pertaining to the impact of X- 
rays on the 21-cm PS are discussed in §(6] In §[7]we consider 
the sensitivity of low frequency arrays to X-ray signatures 
in the 21-cm PS. We summarize our results in §[§1 

Throughout this paper we adopt a concordance cos- 
mology for a flat ACDM Universe, (Sim, Oa, f^b, /i, cs, n) = 
(0.27,0.73,0.046 , 0.7,0.8, 1), con sistent with the WMAP 
constraints from iKomatsu et al.l (200^) . All distances are in 
comoving units unless stated otherwise. 



2 ESTIMATES OF THE POSSIBLE 
CONTRIBUTION OF X-RAYS TO 
REIONIZATION 

In this section we discuss estimates of the X-ray background 
at high redshift. We begin by estimating the X-ray photon 
MFP, and then discuss the X-ray emissivity. High redshift 
star-formation scenarios might lead to two sources of X-rays, 
which we consider in turn. In particular, we estimate the X- 
ray flux from these sources relative to UV flux. Within the 
context of this paper, the role of each of these estimates is 
to motivate different input model parameters for the MFP 
and fractional X-ray contribution to reionization. 

2.1 The X-ray photon mean free path 

As mentioned in the introduction, the long X-ray MFP im- 
plies that they reionize the IGM diffusely via a background, 
rather than locally as for UV photons. The most critical pa- 
rameter in studies of the X-ray contribution is, therefore, 
the MFP. For an X-ray photon of energy Ex (measured in 
eV) the MFP in a n IGM of mean neutra l fraction xni is 
iMadau et al.ll2004l : Furlanetto et al.ll2006h 



: 4:.9x 



^-1/3 



Mpc. 



(1) 



15 J VSOOeV, 

The energy dependence in equation ([TJ means that in order 
to estimate the MFP of X-rays, we must specify their spec- 
tral energy distribution. To estimate the range of possible 
X-ray contributions to reionization we assume a spectrum of 
non-thermal origin with a constant flux per logarithmic in- 
terval [vFi, ~ constant], such that oc 1/-Bx, which yields 

(A) = [ '"^\xF,dEx/ [ F,dE 



4.9xttt 



Mpc X 



dEx 



15 

(gx/300eV)^ 
Ex 



dEx—. (2) 



If we were to assume a ste eper energy spectrum (as in 
iPritchard fc Furlanettd |2007| . for example) , the estimated 
MFP would be smaller. Of course the estimate also depends 
on the range of X-ray energies considered. It is important 
to note that the MFP is necessarily bounded above by the 
horizon scale, c/H, and all MFPs adopted in this paper ad- 
here to this constraint. As an example, we can consider an 
energy range 0.25 keV < -Ex ^ 1 keV. At z = (7,9), and 
assuming aijji ~ 1, this yields (A) « (153, 98) Mpc. 



2.2 Estimates of the X-ray emissivity during 

reionization based on the properties of local 
galaxies 

We next summarize estimates of the X-ray emissivity aris- 
ing from a range of star-formation related sources (including 
supernovae remnants and X-ray binaries) based on observa- 
tions of local galaxies. Assuming that the correlation be- 
tween local star-formation rate (SFR) and X-ray luminosity 
can be extrapolated to high re dshift, the X-ray luminosity 
during reionization is given by l|Furlanettoll2006l ): 



Lx = 3.4 X lO^Vx 



SFR 



IM 



©yr 



erg s 



(3) 



where fx is a normalisation factor that parametrizes igno- 
rance of star formation at high redshift. In this expression 
the SFR can be approximated by 



SFR ^ /, 



1/2 



da 



,1/2 



-nH(2:)mH kgs 



(4) 



where a = {1 + z)~^ is the cosmic scale factor, nu{z) is 
the number density of hydrogen at redshift z, mu is the 
atomic mass of hydrogen, /, is the star-formation efficiency, 
and we have ignored the effect of the cosmological constant 
(appropriate at high redshift). Combining equations @ and 
Q, multiplying by the Hubble time and dividing by nii{z), 
we obtain the luminosity per baryon emitted in X-rays 



ex — tH{z)Lx/nH ■ 



(5) 



Finally, the mean number of X-rays per baryon TVx can be 
found by averaging the number of X-rays per baryon over 
all X-ray energies, which yields 

dex 1 



iVx = 



du- 



dp Ex 



f , dex 1 
= '^J'^liiTE^^ 

= 5.58 X 10"^ afx 



, dex 
dv— — 
dv 

f, dFcou/da I 1 
0.1 0.01 V^x 



(6) 



The factor describing the mean of the inverse energy of X- 
rays is evaluated using the flux density F^, where again we 
assume F^v ^ constant. 



1 

eI 



— F,AEx 
Ex I 



1/ -Bmin — 1/-En- 



F^dEx 



In (i5niax/-Emin) ^ ^ 

For an energy range 0.25 keV ^ Ex ^ 1 keV, this expression 
yields (1/Sx> ~ 1/460 eV. 

It can be shown that for a weakly-ionized IGM only 
about 1/3 of the X-ray energy contributes to ionizatio n 
jShull fc van SteenberdI 19851 : IChen fc Kamionkowskilliooj ) , 
but that e ach X-ray photon pro duces approxima tely 12 
ionizations (|Shull fc van Steenbera .1985. ; Venkates an et al.l 
l200lll . mostly arising from secondary ionizations by pho- 
toelectrons from ionized helium. As reionization proceeds, 
however, the proportion of each X-ray's energy that is de- 
posited into the IGM as heat, increases. In particular, for 
^Hii ^ 1 0%, ionization by second ary electrons becomes in- 
efficient l|Ricotti fc Ostrikerl I2OO4I ) . as the electron-electron 
scattering cross-section becomes large. With these consid- 
erations in mind we estimate the number of ionizing X-ray 
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photons per baryon to be 



iVx,io 



= 12iVx 
« 0.015 /x 



0.1 



,ii/d(lna) ^460eV 



0.05 



Ex 



(8) 



We note that as reionization proceeds, equation ^ increas- 
ingly overestimates the X-ray contribution to reionization. 

With regard to the role of X-rays in the reionization 
of the IGM, it is instructive to compare the number of X- 
ray to UV ionizations. To achieve this we first evaluate the 
star-formation rate of UV sources (which is the same as for 
X-rays here since we are only considering X-rays from stellar 
emission in this subsection) using equation ^ and hence the 
number of ionizing photons per baryon 



A*'uv,io 



0.6 



/esc /* dFcoii/d(lna) 



0.03 0.1 



0.05 



(9) 



where /esc is the fraction of ionizing photons that escape 
from star-forming galaxies. In evaluating equation (|9]) we 
have assumed that 4000 photons arc produced per baryon 
incorporated into stars (Barkana_&_Locb 2001). Using equa- 
tions ((Sj and ([9]) we obtain the relative number of X-ray to 
UV ionizations, which in turn yields the fractional contribu- 
tion to reionization by X-rays 



#^^ 0.025/x^^(^ 



(10) 



As mentioned above, the normalisation factor /x 
parametrizes the uncertainty in the relation between 
the spectral energy distribution produced by low and high 
redshift star formation. Equation (|10p suggests that X-rays 
would need to comprise a large fraction of the energy budget 
when compared with low redshift star formation {fx 3> 1) 
in order for X-rays following directly from star formation to 
make a substantial contribution to the reionization of the 
IGM. 



2.3 X-ray background generated by black hole 
accretion 

An ahernative means for estimating the X-ray flux 
is to consider eflicient accreti on onto pop-Ill r e mnan t 
BHs at high redshift (e.g. iRicotti fc Ostrikej |2004| ). 
In terms of X-ray production, accretion is poten- 
tially more eflicient than star formation. Indeed, sev- 
eral authors have argued that the number of primary 
and secon dary X-ray ionizations per accreted baryon is 
A^x ~ 10*^ (|Ve nkatesa n et al.|[200 ll: 'Ricott i fc Ostrikeri[200i : 



iRicotti et £01 2005: Sh uU fc Ven katcsan 20081). Once again 
we assume that each X-ray photon produces ~ 12 ioniza- 
tions. As a result, the fraction of baryons accreted, cjbh, need 
only be very small to achieve of order one X-ray ionization 
per baryon in the IGM. F or the purpose of illustra tion, we 
follow previous estimates l|Ricotti fc Ostrikerl bOO'j ) by as- 
suming lubh ~ 10~^ — 10~^, with approximately 1/3 of the 
total X-ray energy contributing to reionization. The number 
of X-ray ionizations per b aryon in the IGM then becomes 
(|Shull fc Venkatesanll2008h 



Ax,io 



12 



2.7 X 10-9 erg \ Ex 



460 eV 



(0.01 — 0.1 ions/baryon) 



460 eV ^ 

Ex , 



(11) 



where 0.1 in the numerator is the accretion efficiency and 
the escape fraction and clumping factor are assumed to be 
unity. Using Equations © and (|lip . and assuming that UV 
photon production is dominated by stellar sources, we can 
now estimate the contribution to reionization of X-rays pro- 
duced by black hole accretion relative to UV photons 



Ax,io 



A,, 



(0.02-0.2] 



0.03 



460 eV 



E^ 



/csc.UV 



0.1 



dFcoii/d(lna) 



0.05 



.(12) 



Equation (|12|l demonstrates that the fractional contri- 
bution of X-rays to the ionization of the IGM could be as 
high as 10s of percent during the early stages of reionization. 
This estimate of the number of X-ray ionizations per baryon 
due to accretion is in excess of that expected directly from 
star formation, and demonstrates the plausibility of a non- 
negligible X-ray contribution to reionization. As mentioned 
earlier, one caveat is that theoretical calculations show that 
the free electrons produced by reionization limit the fraction 
of the IGM that can be ionized by X-rays alone to < 20% 
(|Ricotti fc Ostriked[200i ). while observations of the soft X- 
ray backgr ound limit the X-ray contribution to reionization 
to < 50% (|Diikstra et al.l 1200^ 1. It should be noted that in 
this paper the X-ray fraction refers to the fraction of IGM 
ionized by X-rays prior to a particular redshift. In principle, 
this fraction could therefore be as high as 50% until midway 
through reionization. 

Based on the considerations summarized above, we are 
able to identify a plausible range for the contribution of X- 
rays to the reionization of hydrogen. In the modelling pre- 
sented in the remainder of this paper we therefore analyse 
scenarios for which there is a fractional contribution to the 
reionization of hydrogen at redshift z of between 0% and 
50% due to X-rays with mean energies between of 300eV 
and 460eV. These energies yield MFPs that straddle the 
range of values for the characteristic bubble scale during 
the various stages of the reionization era. 



3 SUMMARY OF THE MODEL FOR 
REIONIZATION 

We model fluctuations in the 21-cm signal from neutral hy- 
drogen during reionization using a scmi-rmmerical scheme 
following the work of G eil fc Wvithc ( 2008), which is based 
on the analytic model of lWvithe fc Loebl (|20o3). The an- 
alytic model encodes the preference for galaxy formation 
in overdense regions ||Mo fc Whitelll996l ). in concert with 
the suppression of low mass galaxies in ionized regions re- 
sulting from ionization feedback. The robust prediction of 
this model is that overdense regions reionize first. Individ- 
ual Hll bubbles surrounding clustered sources then perco- 
late into less overdense regions, resulting in the eventual 
overlap of ionized regions throughout the IGM. This be- 
haviour mimi cs that seen in soph isticated numerical sim- 
ulations (e.g. iMcQuinn et al.ll2007l ). In this paper we con- 
sider a model in which the mean IGM i s reionized at z: ^ 6 
l|Fan et al.l |2006| ; iGnedin fc Faiil l2006l : IWhite et al.l |2003| ). 



© 2008 RAS, MNRAS OOO.ITIfTSl 



X-rays & the 21-cm power spectrum 5 



The model differentiates between star formation in neu- 
tral regions, which proceeds in halos above the hydrogen 
cooling threshold Tmin ~ 10* K, and star formation in ion- 
ized regions, which is suppressed below a higher threshold 
Tlon ~ 10^ K due to radiative feedback. The model predicts 
the time dependent ionization fraction by mass Qs^r of a 
particular region on a scale R with overdensity Sr. On scales 
much greater than the ionising photon MFP, this ionization 
fracti on evolves according t o the following differential equa- 
tion (|Wvithe fc Loebll2007l ) 



dt 



iVic 



0.76 



Qs,R. 



dF^M{SR,R,z,Mi 

o 



+ (1 - QS,R.) 



dt 

dFcoll{5R,R,Z,Mr, 



-aBCrin 



1 + 5r 



dt 
D{z) 



-D(Zobs) 



{l + zyQs,R, (13) 



where Mon is the number of ionizing photons per baryon 
in galaxies, ob is the case-B recombination coefficient, C is 
the clumping factor (taken to be 2 in the work presented 
here) and -D(2obs) is the growth factor between redshift z 
and the present. The halo masses Mmin and Mion corre- 
spond to virial masses of temperature Tmin and Tion respec- 
tively. The last term in equation (I13|l describes the recom- 
bination rate of ionized gas. Ionizing photons are assumed 
to be produced at a rate proportional to the collapsed frac- 
tion of mass -Fcoii in halos above the minimum thresholds 
in neutral (Afmin) and ionized (Mion) regions. In a region 
of comoving radius R and mean overdensity Sr{z), the rele- 
vant collap sed fraction is obtained f rom the exte nded Press- 
Schechter (|Press fc Schechteijll974 ) formalism l|Bond et al.1 
M) 



Fcou{Sr, R, z) = erfc 



Sc - Sr{z) 



2 

gal 



(14) 



From the ionization fraction Qs,r we can calculate the cor- 
responding brightness temperature T using 



T{5r,R) = 22 mK 



1 + z 
7.5 



1/2 



{l~Qs.R)il+SR) 



(15) 



4 EXTENSION OF ANALYTIC MODEL TO 
INCLUDE X-RAYS 

We extend the model described in §[3] to account for the 
contribution of X-rays to the reionization of the IGM. The 
X-ray contribution is modelled as six coupled differential 
equations describing the contributions of X-rays and UV 
photons on scales R {Qx,s,r and Quv,s,r respectively) and 
W {Qx,s,\ and Quv,s,\ respectively) to the evolution of the 
ionized fraction on a scale R (Q5,r). Separate rates must 
be computed because at any point in the IGM i) X-rays 
and UV photons cause different levels of radiative feedback, 
ii) they are generated from regions of IGM with different 
scales and overdensities, and iii) the relative contribution 
to reionization from UV and X-ray photons is sensitive to 
radiative feedback on different scales. The equations are: 



dt 



dQuv,s,R ^ dQx,6,\ 



dt 



dt 



(16) 



dQs,\ dQuv,s,x , dQx,s,\ 



dt 



dQuv,s.R 



dt 



+ 



dt 



(17) 



dt 



(1 -Xfrac)Tr^ X 



7i(Qcrit — Qx,5,r) 



Qs,R 



0.76 

dF^aniSR,R,z,Mi 

o 



dt 



+ (1-Q«,fl) 

+Ti.{Qx,S,R — Qcrit) 
(3uV,(5,R 



dFcoll{SR, R, 2, Mmin) 



dt 



dFcon{SR,R,z,Mio„ 



dt 



Qs,R 



(18) 



dQuV,6.\ V- \Nion 

= U — Afrac ) X 



dt 



H(Qcri 



t — '^X,S.X) 



Qs.x 



0.76 

on / 



dt 



-'H{Qx,6,X — Qcrit) 



dFcoll((5(A> , W,Z, Mmin) 



dt 



df;oIl((5(A), (A),Z,Mion) 



dt 



+Bs. 



dQx,s,R 
dt 



iuv,s,x 



(19) 



„ Mon f 1+Sx , 



Ti-iQciit — Qx,s.r) 



0.76 Vl + fe, 

on J 



dt 



-{1-Qs,r) 



~'H{Qx,S,R — Qcrit) 



dF,M{5(x),W,z,Mn 



dt 



dFcoll{S(x), (A), 2, Mion) 



dt 



+Bs,R- 



Qs,R 

dQx,s,x _ 
di ^ 



(20) 



Afrac , X 



0.76 



'H(Qcrit — Qx,S,x) 



rffcoll(<?(A), (A),2,Afion) 



dt 



+ (l-Qi,A) 

—'H{Qx,S,X — Qcrit) 
^X,S,X 



d-Fcoll,x(5(A), (A), 2, Mmin) 



dt 



rf-Fcoll(^(A), (A),Z, Mion) 



dt 



+Bs. 



Qs,. 



(21) 



where the recombination rates on scales R and A, Bs^r and 
Bs.x are given by 



BsM = — QbCtih 
and 



1+5r 



D{z) 



D{Zohs) 



{l + zYQs.. 



(22) 
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D(z) 



{l + zYQs,x. 



(23) 



where Xfrac is the fraction of the total number of ionizations 
that are due to X-rays as estimated in S I2.2I Throughout this 
paper Mon is chosen such that the mean ionization fraction 
at z = 7 is 0.7 and at 2: = 9 is 0.3, which, in the no-X- 
ray case, corresponds to complete reionization at z = 6. 7i 
is the Heaviside step function, which is equal to one when 
the argument is positive definite, and zero otherwise. Fcoii 
is given by equation (|14|l . The term in equation (|20|l that 
is proportional to the ratio of densities on scales A and R 
describes the ionization of a region of overdensity S-r. due to 
X-ray photons generated in a region of overdensity 5\. 

In this paper we consider, for simplicity, only the 
mean scale (A), rather than the full distribution of 
MFPs for X-rays of different energies. Evaluation of 
Jcoii('5(A)5 (-^), z, Mion) requires calculation of the correlation 
between Sa and S(^x) 1 which are concentric. This correlation 
is implicit in the application of equation (|14|l to the semi- 
numerical model in §[5] Ifowever, in our analytic calculations 
we make the further simplifying assumption that the fraction 
of matter that has collapsed into halos and that contributes 
to X-ray ionization in a region of radius R depends on the 
overdensity as 



(FcoI\{Sr,R,z) 
[R'/{\f)F,^ii[5R,R,z) + 
[l-R'/{\f)F,M5o,oo,z) 



R> (A) 



R< (A) 



(24) 



where (5o = is the mean overdensity of the Universe. Equa- 
tion (|24p describes the difference imprinted on large versus 
small scale fluctuations by the addition of X-rays. Firstly, 
if (A) < R then all photons, X-ray and UV, are generated 
inside the region where the fluctuation is measured (i.e. on 
a scale R). When (A) > R, UV photons are generated in- 
side the region R, but X-rays come from both inside R and 
from the region bounded by R and (A). In the latter case. 
X-rays produced inside R are subject to the same bias as are 
UV photons; however X-rays from beyond R are produced 
in a region of different overdensity. In our analytic calcula- 
tions we assume for simplicity that X-rays from beyond R 
are generated at mean density (as already mentioned, this 
assumption is relaxed in §(5]). 

As discussed in S I2.2I approximately 2/3 of the total en- 
ergy emitted in X-rays contributes to heating the IGM, and 
the remaining 1/3 contributes to ionization. As the IGM 
temperature increases, the formation of smaller mass galax- 
ies is suppressed, which should be reflected in the minimum 
scale on which the density field is filtered. We account for 
the effects of X-ray heating by monitoring the tempera- 
ture change of the IGM that is due to X-rays, which we 
assume to be proportional to the ionization fraction due to 
X-rays. When the X-ray ionization fraction reaches the heat- 
ing threshold, Qcrit, defined by 

,TlGM 7? 

104 2 ' 



Qcot = 0.03- 



(25) 



we increase the minimum scale at which the density is fil- 
tered from Tmin to Tion. TiGM ~ 10'* K is taken as the char- 
acteristic IGM temperature above which galaxy formation is 
suppressed. We take Qcrit = 0.03 to be the fiducial value and 
use rj to parametrize our ignorance of the X-ray spectrum, 
as well as the emissivity history. We take Tion = 10^ K to be 



the fiducial value but allow it to vary given the uncertainty 
in Tion when TIgm = 10^ K. This physics is included in our 
model via the Heaviside step functions in equations I18I - I21I 
The simple analytic model described above can be used 
to qualitatively demonstrate the impact of X-rays on the 
statistics of brightness fluctuations during reionization. We 
describe these results in the following subsection before mov- 
ing on to more realistic calculations based on our semi- 
numerical model. We emphasise that throughout this paper 
we assume that the IGM is considerably hotter than the 
CMB, and hence do not consider temperature fluctuations 
as additional sources of 21-cm fluctuations. 



4.1 Results - analytic 

We express results from our analytic calculation via the vari- 
ance of 21-cm brightness temperature within regions of phys- 
ical size corresponding to an observed angle 9. This is given 
by 



2 

0"21 



where 



((T- 



= f dS{T{6,R) 

7% J 



(T)fe 



-S-'/2a 



(T) 



\/2^ 



dST{5, R)e 



-S^ /2a 



(26) 



(27) 



Figures [T] and show a^i (6*) as a function of angular scale 
at 2 = 7 and 2 = 9 respectively. In this example, the mean 
X-ray energy is taken to be 300 eV and fractional X-ray con- 
tributions are shown at the 0%, 10%, 30% and 50% level. We 
reiterate that these contributions reflect the fraction of the 
aggregate ionizing radiation (UV and X-ray photons) that 
are X-rays. 

We first consider the case in which X-ray heating has a 
negligible effect on galaxy formation (Qcrit > 1). The results 
for this fiducial case are shown in the upper left panels of 
Figures [T] and [51 where it can be seen that the variance of 
the 21-cm signal is insensitive to the X-ray contribution for 
all angular scales 9 > 9(^x^ (where 9(^x^ is the angular scale 
corresponding to the MFP). On scales 9 < O/^x) the X-ray 
contribution alters the variance, with the sign of the effect 
depending on the epoch. At 2 = 7, (for the chosen model) 
the variance on scales less than 9(^x) is reduced relative to 
the no- X-ray scenario. This is because the presence of X- 
rays lessens the bias of 21-cm emission toward underdense 
regions (which results from reionization by UV photons that 
are restricted by galaxy bias and short MFP to overdense 
regions of the IGM). It should be noted that the analytic 
calculations cannot describe the effect of ionized bubbles on 
o'2i{9), which imprint a characteristic, redshift-dependent 
scale on the PS. At 2 = 9 the effect of X-rays on a'2i{9) is 
reversed: fluctuations on scales below Oi^x) are enhanced by 
the presence of X-rays. 

These trends can be explained by noting that when the 
average ionization fraction in the IGM is (Q) < 0.5, over- 
dense regions have a h igher brightness temperature than 
underdense regions ((Wvithe fc Morales . 2007 ). For our sim- 
ulations, Q in an overdense region is smaller when X-rays are 
present because the X-ray contribution is proportional to the 
mean overdensity, and hence is reduced relative to the UV 
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Figure 1. The variance of the 21-cm temperature brightness as a function of angular scale at 2 = 7. The X-ray energy is taken to be 
300 eV and X-ray contributions are considered at the 0% (solid), 10% (dotted), 30% (dash) and 50% (dash dot) level. The vertical lines 
(grey) mark the angular scale of the MFP of 300 eV X-rays, which is 7 arc-minutes. 



contribution in an overdense region. Since brightness tem- 
perature increases as Q decreases, we expect to see a relative 
enhancement in brightness temperature in overdense regions 
on scales where the X-ray contribution to reionization differs 
from the UV contribution (i.e. 9 < d{x))- Similarly, in under- 
dense regions Q is greater when X-rays are present and hence 
the brightness temperature is smaller. The combined effect 
of these two factors is to enhance fluctuations in the temper- 
ature brightness if (Q) < 0.5. When (Q) > 0.5 the converse 
holds: on small scales, 8 < 9^x), and relative to the mean 
temperature brightness of the IGM, X-rays suppress fluc- 
tuations by boosting the ionization fraction in underdense 
regions and reducing it in overdense regions. The proper- 
ties of the analytic model when X-ray heating is ignored do 
not vary qualitatively with X-ray energy or fractional X-ray 
contribution. 

The inclusion of X-ray heating of the IGM complicates 
the interpretation of the resulting 21-cm variance. Most no- 
tably, there is no longer quantitative agreement between 
the variance arising from different X-ray fractions on scales 
greater than 9(^x) ■ X-ray heating suppresses the formation of 
small mass galaxies, which increases the clustering bias of 



emission sources. The greater the bias of sources, the earlier 
one expects overdense regions to be ionized relative to un- 
derdense regions. There are, therefore, two competing effects 
due to X-rays: the tendency for X-rays to shift ionizations 
from overdense to underdense regions due to their long MFP, 
and the increase in source bias due to the suppression of low 
mass galaxies. The upper right and lower panels of figures [1] 
and [5] explore the interplay of these two effects by varying 
Qcrit and Tion. In the lower left panels, Qcrit is decreased to 
0.03, such that regions of the IGM that have been heated by 
X-rays to at least 10'* K are filtered on the same mass scale 
as fully ionized regions. The most notable result of reduc- 
ing Qcrit at 2: = 7 is that the suppression of fiuctuations on 
scales smaller than 9^x) is not as prominent as when X-ray 
heating is ignored. At z = 9 this effect is more dramatic, as 
X-ray heating acts to suppress, rather than enhance 21-cm 
fiuctuations. In addition, the results in the panels of figures 
[T]and[2]in which heating has been considered expose a trend 
for fluctuations to increase overall as the X-ray fraction is 
increased. We owe this effect to the precocious ionization 
of overdense regions, and hence of a greater fraction of the 
total hydrogen content of the IGM, due to an increase in 
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Figure 2. The variance of the 21-cm temperature brightness as a function of angular scale at z = 9 (right). The X-ray energy is taken 
to be 300 eV and X-ray contributions are considered at the 0% (solid), 10% (dotted), 30% (dash), and 50% (dash dot) level. The vertical 
lines (grey) mark the angular scale of the MFP of 300 eV X-rays, which is 4 arc-minutes at z = 9. 



the bias of emission sources. As expected, the eflect is more 
pronounced at high redshift compared to epochs closer to 
overlap. 

If we decrease Tion to 3 x 10'' K (lower right panel of 
figures [T] and [2| , with Qcrit = 0.03, the impact of including 
X-ray heating is reduced overall, as the filtering of ionized 
regions is less distinct from that of neutral regions. Finally, 
if we increase Qcrit to 0.1 (upper right panels of figures [1] 
and [2|, the impact of X-ray heating varies depending on 
the fractional X-ray contribution to reionization. For small 
X-ray contributions (~ 10%), the shift in bias to overdense 
regions is the dominant effect of X-ray heating. As the X-ray 
contribution increases, however, the suppression of fluctua- 
tions on scales smaller than 9^x) becomes more significant. 



5 INCLUSION OF X-RAYS IN A 
SEMI-NUMERICAL SCHEME 

Several recent studies of the 21-cm signal during the 
EoR have employed semi-numerical techniques to reduce 
the computational cost of large-scale numerical simu- 



lations (iMesiriger fc Furlanettol l2007l: Zahn et al.l l2007al : 
iGeil fc Wvithd l2008l : IChoudhurv et allboOSh . The advan- 
tage of semi-numerical simulations, as compared to analytic 
models, is the ability to trace the formation and overlap of 
ionized bubbles, which feature prominently in the PS. On 
the other hand, semi-numerical techniques do not capture 
variations in the actual halo count s fro m linear theory. The 
modelling of lMcQuinn et al.l (j2007l ) and lZahn et~all l|2007bf l 
shows this to be a significant effect. 

The semi- numerical model in this paper is based on the 
method implemented bv iGeil fc Wvithd l|2008h . It employs 
a three-dimensional realisation of a Gaussian random field, 
representing fluctuations in the matter overdensity. The sim- 
ulation is performed in a 200^ box, with side length equal 
to 400 and 1800 cMpc for X-rays with mean energies of 300 
and 460 eV (for a flat X-r ay spectrum) respectiv ely, with 
an underlying ACDM PS (jEisenstein fc Hul llQQ^ ) evolved 
linearly to the specified redshift z. The density field is then 
smoothed repeatedly on scales ranging from the pixel size 
to the size of the simulation box. At each scale we calculate 
the ionized hydrogen fraction Q according to equation (|16p . 
The calculation of Q differs from the analytic case in that 
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z = 9, no X-rays 



z = 9, 30% X-rays ^ 
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Figure 3. Slices from 3D ionization maps at redshifts 9 {top) and 7 (bottom) for the fiducial case with no X-ray contribution (left) and 
30% X-ray contribution (right). The colour scale ranges from neutral (black) to ionized (white). The circles in the bottom right corner 
represent the scale of the assumed MFP for 460 eV X-rays. The simulation slices are 1800 cMpc on a side. 



the overdensity on the scale of (A) is calculated explicitly, 
rather than assuming it to be equal to the mean overdensity. 
For this reason, there is no longer a need to differentiate be- 
tween smoothing scales greater and less than (A) . An ionized 
spherical bubble of radius R forms around voxels in which 
Qs,R 1, producing a two-phase (Q = or 1) ionization 
map of the IGM. Voxels that do not form part of an ion- 
ized bubble after the field has been smoothed on all scales, 
are then assigned a partial ionization fraction correspond- 
ing to the value of Q obtained by smoothing the density 
field on the minimum (pixel) scale. Of importance to this 
study, which requires very large box sizes, is that the use 
of equations (I16I)-(I23I) rather than a barrier formalism (e.g. 
iMesinger fc Furlanettol |2007| ) means that partially ionized 
pixels are treated correctly with respect to the ionization 
contribution of unresolved H II regions. 

It should be noted that the ionization field is simulated 
at a fixed proper time, and as a result the corresponding 21- 
cm PS is computed assuming a non-evolving IGM. However, 
in practice, a single observation measures the redshifted 21- 
cm signal over a range of cosmic time corresponding to the 
particular frequency band of the observation. The conse- 
quence of this is that the observed spherically averaged 21- 



cm PS is modified on the largest scales, because the ioniza- 
tion state of the IGM evolves significantly during the light 
travel time corresponding to the frequency band. We do not 
address this evolution in the present study. 

Physically, two principle features are captured in the 
semi-numerical scheme that are not present in our semi- 
analytic formulation. Firstly, the presence of ionized regions, 
which imprint a scale dependent feature on the PS. And sec- 
ondly, the density dependent calculation of X-ray contribu- 
tion. 



6 RESULTS - SEMI-NUMERICAL 

In order to examine the combined impact of H II regions and 
X-ray ionization on the statistics of 21-cm fiuctuations, we 
now consider the results of the semi-numerical realisation of 
the model outlined in §[3l We explore the parameter space of 
average X-ray MFP (A) and the fraction of ionizations due 
to X-rays (refiecting the uncertainty in /x, (A) and /esc) at 
z = 7 and 9. The chosen redshifts straddle the epoch when 
the mass- weighted, average ionization fraction ((3)mass ~ 0.5 
for this model. We consider scenarios in which X-ray heat- 
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ing has a negligible eflFect on galaxy formation, as well as an 
example in which X-ray heating is parametrized as in S I4.1I 
with g^rit = 0.03 and Tion = 3 x 10* at 2 = 7. Table^details 
(Q)mass for each simulation box used in this paper for which 
X-ray heating is ignored. We note that ((5)mass remains con- 
stant at fixed epoch to within a few percent, indicating that 
the scheme conserves photons in the simulations describeclf|. 
Figure [3] shows slices from the three-dimensional simulation 
cubes at z = 7 and z = 9 for the fiducial (no X-ray, no heat- 
ing) case and for an X-ray contribution of 30% (no heat- 
ing). The corresponding mass- weighted ionization fractions 
are 0.77 and 0.32 for the fiducial case, and 0.74 and 0.32 
for the 30% X-ray case at 2: = 7 and « = 9 respectively. 
Also represented as circles on these plots are the size of the 
MFP (« 153 and 98cMpc at z = 7 and z = 9 respectively, 
corresponding to X-rays with mean energy 460 eV) . These 
maps illustrate the primary effect of adding a diffuse X-ray 
background. In particular the maps ai z = 7 illustrate that 
replacement of UV flux in the model with X-ray flux leads 
to smaller H 11 regions, and lower overall contrast between 
21-cm emission from different regions of the IGM. 

It is important to note that the simplifying choice to 
consider only a single X-ray energy, rather than a poly- 
chromatic spectrum, precludes a thorough treatment of how 
cos mological redshift af f ects t he spectral shape. In partic- 
ular iRicotti fc Ostrikej (|2004l ) assert that redshifted hard 
X-rays constitute a large fraction of the ionizing radiation 
duri ng reionization, hen ce affecting the topology of reioniza- 
tion (iRicotti et al.ll2005l ). Moreover, since hard photons have 
a very long MFP, they can travel far from the overdensity 
where they are produced before redshifting to a wavelength 
where they are absorbed. Thus the effect of cosmological red- 
shift is to decorrelate the X-ray ionization with the density 
field even more than expected from the instantaneous X-ray 
SED. These issues can only be addressed in detail by the 
inclusion of radiative transfer. In lieu of radiative transfer 
we have assumed an effective value of MFP. We note that a 
value in excess of that expected from the spectrum provides 
an approximation for the effects of cosmological redshift. 

The results presented in the following subsections are 
influenced by the interplay between three physical parame- 
ters: (i) the MFP of X-rays, (ii) the characteristic scale of 
H II regions, and (iii) the mean ionized fraction in the simu- 
lation box (Q). We consider first the effects of varying these 
three parameters without including X-ray heating. Heating 
is included in the final subsection. We briefly summarize the 
main effects before discussing the different cases at greater 
length with the aid of several examples. 

Both the MFP and the H II region (bubble) scale gener- 
ate a shoulder in the PS, which results from the transfer of 
power between small and large scales. In the case of the X- 

^ For (Q) < 0.5 fluctuations in the temperature brightness are 
dominated by fluctuations in the matter density field; whereas 
for (Q) > 0.5 fluctuations are dominated by HII regions. These 
results, in conjunction with the output of the analytic model, have 
been used to interpret the output of the semi-numerical scheme. 
In particular, at very high-redshift (z > 13) we find that tr^^ (6») 
and the PS are not effected by the X-ray contribution. Similarly, 
(Q) ^ prior to the formation of bubbles, we find that 
the semi-numerical scheme produces results very similar to the 
analytic model. 



Table 1. Table of mean mass- weighted ionization fraction for 
different simulation parameters when X-ray heating is ignored. 



z 


side length 


(Ex) Rat eV 


(A) cMpc Xfrac 


(Q>mi 


7 


400 


300 


19 


0.0 


0.81 


7 


400 


300 


19 


0.1 


0.80 


7 


400 


300 


19 


0.3 


0.78 


7 


1800 


460 


153 


0.0 


0.70 


7 


1800 


460 


153 


0.1 


0.70 


7 


1800 


460 


153 


0.3 


0.69 


9 


400 


300 


12 


0.0 


0.32 


9 


400 


300 


12 


0.1 


0.32 


9 


400 


300 


12 


0.3 


0.32 


9 


400 


300 


12 


0.5 


0.31 


9 


1800 


460 


98 


0.0 


0.29 


9 


1800 


460 


98 


0.1 


0.29 


9 


1800 


460 


98 


0.3 


0.29 


9 


1800 


460 


98 


0.5 


0.29 



ray MFP, the transfer is across scales comparable with the 
MFP. Conversely, the shoulder generated by H II regions re- 
sults from the transfer of power from small scales to scales 
exceeding the size of ionized bubbles. When these two scales 
are comparable, it is difficult to attribute the features of the 
PS to a single physical process. On the other hand, when the 
MFP and H II region scales differ substantially, the resulting 
PS depends on (Q). The presence of X-rays decorrelates the 
density and brightness temperature, most notably on scales 
smaller than the MFP. Thus, as discussed in § 14.11 when 
(Q) J; (^)0-5, the brightness temperature of overdense re- 
gions is enhanced (suppressed) relative to that of underdense 
regions. In contrast, as discussed in the same section. X-ray 
heating of the IGM increases the bias of sources by suppress- 
ing low mass galaxies, resulting in overdense regions being 
ionized first, even in the presence of a large X-ray fraction. 

6.1 Models without X-ray heating 

6.1.1 Case 1: the mean free path is smaller than the 
bubble scale 

Figure |4] presents results for 21-cm fluctuations ai z — 7 
for (A) — 19cMpc and (A) — 153cMpc (corresponding to 
(Ex) = 300 eV and (Ex) = 460 eV respectively for a flat 
spectrum). In the upper panels of figure |4] we show the 
PS of 21-cm fluctuations computed from the semi-numerical 
model for 300 eV X-ray photons at z = 7. At this late stage 
of reionization the MFP of 300 eV photons is much smaller 
than the H II bubble scale, which means that the X-ray pho- 
tons that are contributing to ionizing the bubbles are being 
produced within the bubbles. For this reason, the bubble 
scale is not affected by X-rays in this illustrated 
by the PS shown. However, we find that on scales smaller 
than the MFP, the presence of X-rays acts to suppress fluc- 
tuations in the brightness temperature. The different lines 
represent X-ray contributions to reionization of between 0% 
and 30% (the case of 50% X-ray contribution is not shown 
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Figure 4. Left: The power spectrum {left) of fluctuations in the 21-cm brightness temperature at 2; = 7 for X-ray contributions of 0%, 
10% and 30%, and mean X-ray energy of 300 eV {top) and 460 eV {bottom). The vertical hne {grey) marks the position of the MFP. 
Right: Corresponding fluctuations relative to the fiducial no X-ray case. Also plotted for comparison is the estimated sensitivity of the 
MWA (upper left panels, solid grey) and MWA5000 (upper left panels, dashed grey), and the sensitivity relative to the PS (upper-right 
panels) assuming 1000 hr integration on 1 field, and bins of width Afe = fc/10. Tion = 10^ K, Qcrit > 1- 
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Figure 6. As per Figure g] with Tion = 3 X 10'' K, Qcrit = 0.03. 



because this would correspond to X-rays reionizing ~ 40% 
of the IGM, which exceeds the hmits discussed in S I2.3|) . Also 
shown (upper-right panel) is the absolute value of the fluc- 
tuation due to X-rays relative to the fiducial case (with zero 
X-ray contribution). 

Figure |4] demonstrates that as the X-ray fraction in- 
creases, the suppression of fluctuations below is more 
pronounced, in agreement with the analytic model. 

6.1.2 Case 2: the mean free path is comparable to the 
bubble scale 

In the upper panels of figure [5] we repeat this analysis at 
2 = 9, when (Q) < 0.5. In this case, the scale of the X- 
ray MFP is comparable to the characteristic bubble scale, 
and hence the prominence of the bubble shoulder is dimin- 
ished in the presence of X-rays. As the relative X-ray con- 
tribution increases, the bubble shoulder is less visible. As 
a result of the suppression of the bubble shoulder, fluctua- 
tions on scales larger than the MFP are reduced. Moreover, 
on scales smaller than the X-ray MFP, fluctuations are also 
suppressed, albeit to a lesser extent. Hence this result con- 
trasts with the z — 7 case in the upper panels of Figure |4j 
in which suppression of fluctuations occurred only on scales 
smaller than the MFP. Once again, we invoke the decorrela- 
tion of density and ionization to explain this result. That is, 
X-rays reduce the biased reionization of overdense regions by 
UV photons, and hence reduce the contrast between ionized 
regions in the IGM and the neutral component of IGM. 
The upper panels of flgure [5] show the case of a MFP 



that is comparable to the bubble scale early in reioniza- 
tion. However this condition is also met for a larger MFP 
later in reionization, as shown in the lower panels of fig- 
ure U where the PS (left) of 21-cm fluctuations at z = 7 
for {Ex) = 460 eV is presented. The different lines represent 
X-ray contributions to reionization of between 0% and 30%. 
Also shown (right-hand panel) is the fluctuation relative to 
the fiducial case. These results show that the maximum sup- 
pression of fluctuations occurs on the scale of the MFP. As 
the X-ray fraction increases, the shoulder due to H II regions 
that is present in the flducial PS (solid) becomes less promi- 
nent. As mentioned previously, this is because this feature 
is generated by the transfer of power from small to large 
scales following the formation of H II regions, an effect which 
is washed out by the diffuse ionization due to X-rays. 

The lower panels of flgure |3] demonstrate that as the X- 
ray fraction increases, the suppression of fluctuations below 
(A) is more pronounced. In accordance with the analytic 
results presented in S I4.1I we flnd that the suppression of 
fluctuations due to X-rays is greatest at a scale slightly less 
than (A). This is most easily seen in the plots of ratio relative 
to the flducial case in the lower right panel of flgure |4] 

6.1.3 Case 3: the mean free path is larger than the bubble 
scale 

In this subsection we explore the impact of increasing the 
X-ray MFP, such that the MFP is much larger than the char- 
acteristic scale of Hll regions. In the lower panels of Figure 
[5] we repeat our analysis at z = 9, when (Q) < 0.5. Most 
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notably, the X-ray contribution magnifies fluctuations in the 
21-cm signal, rather than suppressing them. As discussed in 
§ 14.11 this reversal can be attributed to the contrasting effect 
on the brightness temperature of underdense and overdense 
regions at epochs either side of (Q) ~ 0.5. As at lower red- 
shift, the scale of maximum departure from the fiducial case 
is slightly less than 6{\)- However, unlike the low redshift 
case, the ratio of the power spectra asymptotes to 1 at large 
scales, which we attribute to there being very few ionized 
bubbles at high redshift. 

We find that there is the greatest qualitative agreement 
with the analytic results presented in figure [2] in this case, 
which we attribute to the lesser importance of Hll regions 
on scales much larger than their characteristic size. Con- 
versely, we attribute the discrepancy between the analytic 
and semi-numerical results closer to the end of reionization 
to additional fiuctuations due to finite numbers of H II re- 
gions late in reionization, which are not captured by the 
analytic model. 



6.2 Models with X-ray heating: an example 

One possible parametrization of X-ray heating is to set an 
ionization fraction due to X-rays (Qcrit) above which small 
mass galaxy formation is suppressed. The degree of suppres- 
sion of galaxy formation is governed by the filtering scale of 
ionized regions (Tion). Figure [6] shows the results of simu- 
lations at z = 7 for Tion = 3 x 10'* and Qcrit = 0.03. As 
in earlier figures, we consider 300 eV (upper) and 460 eV 
(lower) X-rays, and fractional X-ray contributions of be- 
tween 0% and 30%. For the no- X-ray case the reduction 
of Tion from 10^ to 3 x 10^ K means that the prominence 
of the bubble shoulder is diminished in comparison to the 
case presented in the upper panels of Figure In addition, 
when the X-ray MFP is smaller than the bubble scale, the 
inclusion of X-rays increases fluctuations on all scales. On 
the other hand, as discussed in S I4.1I when the X-ray MFP 
exceeds the bubble scale (lower panels), the profile of 21- 
cm fluctuations is governed by two competing effects: an 
increase in the clustering bias of emission sources due to the 
suppression of low-mass galaxy formation; and a decrease in 
the bias due to the relatively long MFP of X-rays shifting 
ionizations from overdense to underdense regions. Which of 
these two effects dominates depends on the fractional X-ray 
contribution to reionization. For example, when the X-ray 
contribution remains below ^ 20%, the increase in fluctua- 
tions due to the increased source bias governs the departure 
from the no-X-ray case. As the X-ray contribution increases, 
however, fluctuations are suppressed either side of the scale 
of the MFP, due to the long X-ray MFP. 



7 SENSITIVITY OF LOW-FREQUENCY 

ARRAYS TO THE EFFECT OF X-RAYS ON 
THE 21-CM PS 

Before concluding, we compute the sensitivity with which 
the effect of ionization by X-rays on the shape of the PS 
could be detected using forthcoming low-frequency arrays. 
To compute the sensitivity AP2i(k) of a radio interferom- 
eter to the 21-cm PS , we follow the procedure ou tlined by 
iMcQuinn et all (120061 ) and lBowman et all l|2006al ') (see also 



IWvithe et al.ll2"008h . The important issues are discussed be- 
low, but the reader is referred to these papers for further 
details. The sensitivity to the PS comprises components due 
to the thermal noise, and sample variance within the finite 
volume of the observations. We consider observational pa- 
rameters corresponding to the design specifications of the 
MWA, and of a hypothetical follow-up to the MWA (termed 
the MWA5000). In particular the MWA is assumed to com- 
prise a phased array of 500 tiles. Each tile contains 16 cross 
dipoles yielding an effective collecting area of = 16(AV4) 
(the area is capped for A > 2.1m). The physical area of a 
tile is ^tiie = 16 m^. The tiles are distributed according to 
a radial antenna density of p(r) oc r~^, within a diame- 
ter of 1.5 km and outside of a flat density core of radius 
18 m. The MWA5000 is assumed to follow the basic design 
of the MWA. The quantitative differences are that the tele- 
scope is assumed to have 5000 tiles within a diameter of 
2 km, with a flat density core of 80 m. In each case we as- 
sume one field is observed for 1000 hr. Following the work of 
McQuinn et al. (2006) we assume that foregrounds can be 
removed over 8 MHz sub-bands, within a 32 MHz bandpass 
[foreground removal therefore imposes a minimum on the ac- 
cessible wavenumber of knVm '~ 0.04[(1 + 2;)/7.5]^^Mpc~^]. 

The sensitivity to the 21-cm PS per mode may be writ- 
ten 



5P2i(k,9) = 



Ti 



D^AD 



Ai/tint n(kj_) 



+ P2i(k,e), 



(28) 



where D is the comoving distance to the centre of the survey 
volume which has a comoving depth AD. Here n(k±) is 
the density of baselines which observe a wave vector with 
transverse component k±, and 9 is the angle between the 
mode k and the line of sight. The thermal noise component 
(first term) is proportional to the sky temperature, where 
Tsky ~ 250 (iy^)^'^K at the frequencies of interest. The 
second term corresponds to sample variance. The overall 
sensitivity is 

AP21 (fc, 9) = SP21 (k, 9)/^N4k,9), (29) 

where Nc(k,9) denotes the number of modes observed in a 
fc-space volume d^fc (only modes whose line-of-sight com- 
ponents fit within the observed bandpass are included). 
In terms of the fc-vector components k and 9, Nc — 
2TTk'^sm9dkd9V/(2TT)^ where V = D'^AD(\'^/A^) is the ob- 
served volume. Taking the spherical average over bins of 9, 
the sensitivity to the 21-cm PS is 



[AP2l(fc)]2 



[AP2l(fe,0)]2 



The spherically averaged sensitivity curves for the 
MWA and MWA5000 (within bins of Afc = fc/10) are plot- 
ted as the solid grey and dashed grey lines in the upper left 
panels of each of Figures r4l6l The sensitivity as a ratio of the 
fiducial PS (AP21/P21) is plotted in the upper right panels 
of these Figures (again as solid grey and dashed grey lines) . 

Comparison of the sensitivity with the fiuctuations in- 
duced in the PS by X-ray ionization will be measurable by 
the MWA in the cases of long MFP. This is because in these 
cases the fluctuations are largest on scales comparable to 
the MFP, while the low resolution of the MWA results in 
its sensitivity being best on large scales (~ O.lMpc"^). In 
the examples of smaller MFPs shown, the MWA will not be 



© 2008 RAS, MNRAS OOO.mfTSi 



14 Warszawski, Geil & Wyithe 



sensitive to fluctuations due to X-rays, either because fluctu- 
ations are negligible at scales near ~ 0.1 Mpc~^ (figure[4]) or 
because the fiducial PS is barely detected and the presence 
of X-rays suppresses, rather than enhances the fluctuation 
amplitude (flgure[5}. In these cases, comparison with the 
sensitivity curve of the MWA5000 shows that detection of 
the fluctuations in these cases will be possible using an array 
with larger collecting area. 

We note that for the cases with large MFPs, the PS in 
this paper have been plotted out to very large scales. How- 
ever, in practice the requirement that bright foregrounds be 
removed from the 21-cm observations prior to construction 
of the PS is likely to restrict observations to wave numbers 
k > 0.03 - 0.05 Mpc"^ The lower panels of figures g] and [5] 
therefore illustrate that foreground removal will likely ren- 
der the PS undetectable on scales comparable to and a bit 
larger than the X-ray MFP in these cases. As a result, the 
crossover from scales most effected by X-rays to scales less 
effected by X-rays will not be detected. However, as noted 
above, these figures also show that the enhancement (or sup- 
pression, depending on epoch) of the PS on scales below the 
X-ray MFP, which is sensitive to the fractional X-ray con- 
tribution, will still fall within observable scales. Thus we 
expect that measurements of the contribution of X-rays to 
reionization will be accessible to observations of the 21-cm 
PS for all MFPs, while measurements of the mean X-ray 
energy will only be possible for sufficiently low energies. 



8 SUMMARY 

The presence of an ionizing background of X-rays imprints 
a scale dependent signature on the 21-cm power spectrum 
during reionization. The impact of X-rays on the topology of 
reionization depends on the redshift at which one observes, 
as well as the mean energy of the X-rays and their fractional 
contribution relative to UV emission. The findings of this 
investigation can be summarized as follows: 

• The effect of X-rays on the statistics of 21-cm fluctu- 
ations can be prominent on any scale, rather than only on 
scales comparable to the MFP as might be naively expected. 

• In a mostly neutral IGM, X-ray ionizations enhance 21- 
cm fluctuations. Conversely, in a mostly ionized IGM, X-rays 
suppress fluctuations, provided that X-ray heating does not 
suppress low mass galaxy formation. 

• ionization by X-rays modifies, and in some cases 
smoothes out, the signature of Hll regions in the 21-cm PS. 
If reionization of the IGM included a non-negligible contri- 
bution due to X-rays, then the modification of the PS shape 
introduces a feature at a level that is comparable to that 
introduced by H II regions. 

• The details of PS modification depend on the relative 
scales of the H II regions and MFP, and on the overall ion- 
ization fraction of the IGM. 

• The increase in source bias due to suppression of low 
mass galaxy formation by X-ray heating counterpoints the 
decorrelation of ionizing photons and density due to the long 
X-ray MFP. 

• The MWA will have sufficient sensitivity to detect the 
modification of the PS due to a 10-30% contribution to 
reionization by X-rays, so long as the MFP falls within 



the range of scales over which the array is most sensi- 
tive (~ O.lMpc-i). In cases where the MFP takes a much 
smaller value, an array with larger collecting area would be 
required. 

• In some scenarios, the requirements of foreground re- 
moval will limit the observability of the effect of X-ray ion- 
ization on the PS at scales larger than the MFP. However, 
in all cases, the magnitude of the signature imprinted on the 
21-cm PS at small scales is sensitive to the relative X-ray 
fraction, and will fall within the observable range. 

In conclusion, our study shows that the details of the 
PS shape depend on the typical X-ray photon energy, while 
the magnitude of any modification to the PS is determined 
by the relative contribution of X-rays to reionization. The 
possible contribution of X-rays therefore has the potential to 
substantially complicate analysis of the 21-cm PS. On the 
other hand, precision meeisurements and modelling of the 
21-cm PS promise to provide a method for investigating the 
role and contribution of X-rays during reionization. 
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